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Abstract 

The free energy of the chiral Potts model has been obtained in two ways. The first 
used only the star-triangle relation, symmetries and invariances, and led to a system of 
equations that implicitly define the free energy, and from which the critical behaviour 
can be obtained. The second used the functional relations derived by Bazhanov and 
Stroganov, solving them to obtain the free energy explicitly as a double integral. Here 
we obtain, for the first time, a direct verification that the two results are identical at 
all temperatures. 

1 Introduction 

The chiral Potts model is the first model found that satisfies the star-triangle relation 
but does not have the "difference property". This means that the model and its 
properties cannot be simply expressed in terms of elliptic functions. Its free energy 
was first investigated by showing that it satisfied certain partial differential equations 



*This work was performed while the author was a Visiting Miller Professor at the University of California 
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||] . Later an explicit formula was obtained as a double integral [^, [3| , using the transfer 
matrix functional relations derived by Bazhanov and Stroganov Q and discussed by 
the author || . These two results were shown to be consistent in the scaling region near 
criticality ||, but their equivalence has not been verified in general. Here we do this 
and discuss the properties of the functions x p ,y p that play a key role in |fflj. 

The model is defined in Ref. @. Here we shall consider only the model on the 
square lattice, but because of the star-triangle relation corresponding models can also 
be defined on the triangular and honeycomb lattices. It is shown in |Q] how their free 
energies can easily be obtained from that of the square lattice. 

Here we use the notation of [|J. Let k, k' be two real numbers, between and 1, 
satisfying 

k 2 + k' 2 = l . (1) 
Define variables 9 P , 4> p , u p , v p , related to the a p ,b p ,c p ,d p of jjj] and to one another by 

6 p — 6 ' b p /c p , e ^ p — a p /d p , 

u p = N(9 P + 4> p )/2 , v p = N(6 p - <t> p )/2 (2) 
Then from eqn. (9) of Jjj, u p and v p are related by 

sin tip = k sin Up . (3) 

Thus if the constant k is given, then any one of the parameters 9 P , (j) p , u p , v p 
specifies the other three, to within multiple, but discrete, values. Here we shall regard 
u p as the primary (real) variable, and usually confine it to the interval 

< u p < tt . (4) 
Then there is a unique choice for v p such that it is real, satisfying 



- 7T/2 < V p < 7T/2 . 

The parameters 6 P , 4> p are then uniquely defined by (||). 



(5) 



We also define a function 

T(9,n) = cos' 



N9 



-n/N n 



tt(2j-1) - 
2 2N 



(6) 



Consider a square lattice of L sites. At each site i there is a "spin" <jj which takes 
values 0, 1, . . . , N — 1. Adjacent horizontal sites i,j (with j to the right of i) interact 
with Boltzmann weight function W pq (ai — dj), and adjacent vertical sites i,j (j above 
i) with Boltzmann weight function W pq (ai — Cj). The functions W, W are defined (to 
within normalization factors) by 



W pq (n)/W pq (0) = T{e q -4> p ,n)/T{e p -4> q ,n) , 
W pq (n)/W pq (0) = T(^-</, q + C,n)/T{0 q -0 p -C,n) , (7) 

where 

( = ir/N . (8) 

The above relations ensure that the weight functions are periodic of period N: 
W pq (n + N) = W pq (n), W pq {n + N) = W pq (n). 
Define 




A pq = [det N W pq (i-j)\ (9) 

i.e. is the determinant of the cyclic N by N matrix with entry W pq (i — j) in 

position 

Let Z be the partition function of the model. Then from (2.16), (3.22) and (3.42) 
of [||, the partition function per site is 

Z l l L = Ppq A pq e~ A v« (10) 

so Apq is the dimensionless free energy per site in the normalization in which p pq , A pq 
are both one. This is the function A pq of Q. 

Now we consider the result (21) of 0. There we used the normalization p pq = A pq = 
1, and V {t q ,\ q ) l / L is the partition function per site, so the lhs of (21) is — ANA pq . 

We do need to reconcile the notations. The variables u p , u q of Q are the same as 
those above, except that they lie in the interval — tt < u p , u q < —tt/2. We replace them 
by u p — 7r, u q — 7r (and negate v p , v q ). This is equivalent to replacing the rapidities p, q 
by R~ 1 p, R~ l q, which in turn is equivalent to merely rotating the lattice through 180°, 
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so does not change the free energy, but does ensure that the new u p , u q each lie in the 
interval (m, so we can directly compare the result of [Q] with the equations of 
Making these substitutions, eqn. (21) of |§] gives, for < u p ,u q < ir/2, 

/oo 
e 2 ^ u "- u ^{s(P)[G p (f3)G q (-P) + cosech 2 /?] + 
-co 

t{p)[G P (P) + G q (-0)]}d0 , (11) 

where V denotes the principal-value integral, 

r = (N-l)/(2N) , 

J p = [1 + k 2 + 2k cos(u p + v p )]/k' 2 = sin 2 (u p + v p )/(k' 2 sin 2 u p ) , (12) 
s((3) = [N sinh^ cosh(iV - l)/3 - smhN(3]/[AN(3 sinhiV/3] , (13) 
t{p) = sinh(7V- 1)P/[4P sinhiV/3] , 



and 

q m\ = p / . (14) 

P vr 7-oo sin(u p + ix)(l + k 2 sinh 2 x) 1 / 2 

Note that negating x is equivalent to complex conjugating the integrand in (14), so 



G p {(3) is a real function. Its value at f3 = is 

G p (0) = l-2t> p /vr . (15) 



In fact the integrals are convergent provided < u p ,u q < it, so (|TT|) should hold 
throughout this larger domain, each side being analytic in u p ,u q . For u p < u q the 
Boltzmann weights are real and positive and A pq is the true free energy. For u p > u q 
it is its analytic continuation. By negating (3 it is easily seen that A qp = —A pq , in 
agreement with (3.40) - (3.42) of ffl: this is the "inversion relation". 



2 Derivatives of G p (/3) 

The essential point of the method of [0 is that (because of the star-triangle relation and 
the fact that W pq (n), W pq (n) each depend on only two parameters, rather than three) 
certain derivatives of A pq can be expressed in terms of "single-rapidity" functions that 
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depend on only u p or u g , but not both. (They also depend on k.) We now show that 
this follows directly from (|TTD . 

To differentiate G p {(3) with respect to u p , keeping k fixed, first change the integra- 
tion variable to y, where x = y + iu p . Then differentiate, using (||), and then change 
back to x as the integration variable. The result is that the derivative of G p {(3) is equal 
to the rhs of dl4|), but with an extra factor in the integrand: 

k 2 sin u p cos u p ik 2 sinh x cosh x 
cos 2 v p 1 + k 2 sinh 2 x 

(The first term comes from differentiating cosv p , the second from differentiating the 

integrand.) 

Expanding, this factor becomes 

k 2 sin(u p + ix) [cos u p cosh x + ik' 2 sin u p sinh x] 



cos 2 v p (1 + k 2 sinh x) 



(16) 



We see that sin(it p + ix) appears in the numerator. This cancels the same term in 
the denominator in $L$ ) , leaving 

—G p {(3) = [cos Up A((3) — k sin «p B{fl)\ , (17) 

ou„ cosv„ 



where 



1 f°° cos (2f3x/ir) cosh x dx 



m = l ^wy^+^u (18) 



(1 + k 2 sinh 2 x) 3 / 2 
1 f°° sm(2/3a;/7r) sinh x dx 
(1 + /c 2 sinh 2 x) 3 / 2 

We can also differentiate G p (f3) with respect to fe, keeping u p fixed. This time the 
extra factor in the integrand is 

k sin 2 u p k sinh 2 x k sin(u p + ix) sin(u p — ix) 



cos 2 v p 1 + k 2 sinh 2 x cos 2 v p (1 + /c 2 sinh 2 x) 

Again sm(u p +ix) appears in the numerator, so cancels out of the full integrand, leaving 

— G p (f3) = [sinup A((3)+ cos up B ((3)} . (19) 

Ofv COS Vp 

Note that A(f3), B(f3) are independent of u p . Although the rapidity parameter u p is 
locked into the integral expression (|l~4] ) for G p (j3), it can be taken outside the integrals 
occurring in the derivatives of G p {(3). This is the key to this calculation: indeed, most 
of what we do from now on is merely elementary differentiation and straightforward 
(if cumbersome) algebraic manipulation. 



3 Derivatives of A pq 

First differentiate A pq with respect to k, keeping u p and u q fixed. Noting that 

— log J p = 2 cos Up/ (k' cosvp) , 

OK 

and using (|19D, we obtain 

d , ^ B sin u„ + m„ cos u„ — £ q sin u„ — m„ cos u. 

'J p COS 

where 



. ^p OilJ. U/g |^ / / p *^0/0 OILL Uup II Lq tUO UUp , , 

"^7^p<j = ; j (20) 

ok k cos v v cos ?; 



£ p = k 2 L p cos v p , m p = T(k/k' ) cos i> p — k 2 M p cos v p , (21) 
L p , M p being the integrals 

/oo 
e (3-2(3u p /n [ a (p) Gp {p) + t{(3)] A(J3) dP , 
-oo 

M p = I eP- 2 ^[s([3)G p ([3)+t(l3)}B(l3)d(3 . (22) 



du~ l ° gJp = 



DO 

Secondly, differentiating A pq with respect to u p , u q and summing, noting that 

d , . 2fcsinu„ 
ip 

we obtain 



/ d ^ d \ ^ £ p cos ti g — k' 2 m p sin n g — ^ g cos u p + k' 2 m q sin u p 
I chx p P9 cos Vp cos v q 

The contributions of the term (t/2) log( J q / J p ) have been incorporated into m p . 

Now define functions x p , y p so that 

^ p = (e — x p ) cos Up + k' 2 y p sinu p , m p = (e — x p ) sin u p — y p cos u p . (24) 

Here e is a "constant" - independent of u p but dependent on A;. It is included with x p 
for later convenience. 

Solving these equations for x p , y p gives 

x p = e — {ip cos Up + k' 2 m p sin « p ) / cos 2 v p , 

2/p = (^p sin itp — m,p cos m p ) / cos 2 v p , (25) 

while substituting them into (f20|), (p3|) gives 

d A _ sin(n p + « g )(x g - x p ) + C pg(y g - y p ) 

'Jp COS V q 



A = —™ ' ^ ~" ' ™™ »*> , (26) 
ok k cos v„ cos v„ 
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9 _d_\ A = Cpq(x q - x p ) - k' sm(u p + u q )(y q - y p ) 
dup du q J pq cos v p cos v q 

where 

c w = cos tip cos u q — k sin -u p sm . 

These are precisely the equations (3.45), (3.46) of [jD], so we have verified that the 
double integral expression obtained in [§] has derivatives of this form, and have (for 
the first time) obtained explicit expressions for the single-rapidity functions x p , y p in 

i- 

Note that these functions x p ,y p are not to be confused with the simple ratios 
a p / dp, bp/cp of the original rapidity parameters elsewhere referred to as x p: y p . 

We shall find it helpful to regard G p (f3), . . . , x p , y p as functions G(u p , (3), . . . , 
x(up),y(u p ) of the variable u p , defining cosv p as the positive square root of (1 — 
k 2 sin 2 u p ) for u p real. The functions are defined by the above integrals for < u p < it. 
Outside this interval and in the vicinity of the real axis they are defined by analytic 
continuation. 

From now on we drop the suffix p and regard u as an independent variable; v is a 
dependent variable defined by (^), in particular cos?; = — fc 2 sin 2 u). We refer to 
quantities that are independent of u as "constants" : they may still (and usually do) 
depend on k. 

With these definitions, G(u, (3), . . . , x(u),y(u) are real analytic functions of u for u 
real. In the complex plane they have branch points at 

u = (n - l/2)vr ±cosh _1 (l/fc) (28) 

for all integers n. 

Note that the function G(u,(3) of this paper is not the function G(u) of Q. Nor 
should the function x = x p be confused with the integration variable x in ( |i~4| ) - ( |l9| ) . 
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4 Differential equations for x, y 

An essential step in the working of was the derivation of the pair of coupled partial 
differentiation equations (4.12) for x(u) and y{v): 

2x + I — k— I cos 2 v = a — (a + 5) sin 2 u , (29) 

\ou ok J 

-fe 2 xsin2u + ( — — h kk' — — 2y ) cos 2 w = |(fc' a + 5) sin 2u + rj cos 2 u . 



K du dk 

Here a, 5, rj are unknown constants. (This r\ is the /3 of Q.) 

Much of Q is concerned with evaluating these constants, but we can now obtain 
explicit expressions for them. Using (|2l]) - ( |25| ) to express x, y in terms of L, M, and 
ignoring for the moment the contribution of the constant term e in the definition of 
x(u), we find that the left-hand sides of the two equations ( |29j ) are, respectively, 

k 2 (R cos u + S sin u) cost" , k 2 (k' 2 Rsmu — Scosu) cosv , 



where 



T BL dM 
R = L + k— + — 

ok ou 

S = ^ + (3k 2 -l)M-kk' 2 ^ . (30) 
ou ok 



From the definitions (18) we can establish 



A(/3) + k^jjp--2pB(j3)/* = , 

-2/M(/?)/tt + (3k 2 - l)B{f3) - kk' 2 ^^- = . (31) 

ok 



Substituting the definitions (|22|) of L(u) and M(u) into (|3C|) , and noting that s(/3), 
t(/3) are independent of both u and k, it follows that the only contributions to R and 
S come from the derivatives of G(u,ft). Using ( |i~7D and (|T^), we obtain 

i? = A; 2 Ucosn/cost> , S = k 2 k f2 V sinu / cosv , (32) 

where 

/oo 
s(0)A(l3)B(/3)d/3 . (33) 
-oo 

(There are also terms involving the integrals of s(/?)^4 2 (/?) and s((5)B 2 (f3), but these 
have odd integrands, so vanish.) Hence the left-hand sides of ( p9|) are 

fe 4 U[-l + (2-A: 2 )sin 2 n] , -A: 4 A: ,2 U sin 2u , 



respectively. These are indeed of the same form as the rhs, and at this stage a, S, r\ 
have the values —k^V, —k^k' 2 V, 0, respectively. 

We are still not quite in a position to directly compare our results with There 
it was pointed out that that one can add arbitrary constants to x(u) and y(u) without 



changing (^), (27), and this freedom was used to ensure that t] = and a + 5 = 0. We 
already have rj = 0, which means we do not have to adjust y(u). We do have to adjust 
x{u): we can do this by a suitable choice of e, which parameter gives contributions 



2e — A; cos v de/dk, —kesin2u to the left-hand sides of (p9|), and hence contributions 
2e — k de/dk, — 2e + kk' 2 de/dk, to a, 5 and rj. Altogether we therefore have 

a + 5 = -k 4 V - k 4 k' 2 V - k 3 de/dk , (34) 

so to make a + 5 = we choose 

de/dk = -k(l + k' 2 )V . (35) 



If we define 



/oo 
0' 1 s((3)A 2 (/3)df3 
-oo 

Y = {itk 2 k' A /2) / /T 1 s(p)B 2 ((3)d(3 , (36) 



then from ( pT[ ) we can verify that 



% = kV > ^ = - kk ' 2y • ( 37 ) 
dk dk 

so we can take 

e = Y-X . (38) 

This choice ensures that 

lim x(u) = , (39) 

u— >±ioo 

and the final adjusted values of the constants in ( ]29| ) are 

a = -5 = 2[k 2 k' 2 V + Y - X] , rj = . (40) 



The constant A that figures prominently in [Q] is 

k^d_ (E 
kdk \~k 2 



k' d fk'a\ 

A = -T3lk • ( 41 ) 



5 Properties of the functions 



Periodicity 

We should perform what checks we can to see if the functions x(u),y(u) have the 
properties used in Q. Here we investigate their periodicity relations under the mapping 
u —* u + ir. Here u p and u are the same variable. 

First we have to analytically continue G(u,/3). As u p is increased beyond ir in fll4|), 
the contour of integration has to shifted above the real axis so as to remain above the 
pole at x = i{u p — ir). The effect is the same as leaving the contour on the real axis 
but adding the contribution from a small circle, traversed clockwise, round this pole. 
This contribution is e@ . 

Thus for ir < u < 2ir, G(u, 13) is defined by ([14]), but with an additional term e@ on 
the rhs. It follows that for < u < ir, 

G(u + tt,i3) = -e 213 G(u, 13) + 2 e 13 . (42) 

Similarly writing L p ,M p as L(u p ), M(u p ), and noting that 

a(P) + t(J3) cosh/3 = t/(2(3) , (43) 

it follows from (E22|) that 



L(u) + L(u + tt) = tV e- Wu,n l3~ l A{(3)d(3 , 

J — oo 
/•oo 

M(u) + M{u + tt) = t e~ 2pu/n f3~ 1 B{P)df3 . (44) 



oo 



The factors cos(2(3x/tt), sin(2/3a; / n) in ( |lq ) can be written as the real and imaginary 
parts of exp(2i/3x/7r). Doing this, substituting the resulting expressions for A(/3), B{(3) 
into (|44]) , differentiating with respect to u or integrating by parts in the expression for 
B(f3), we can arrange that the (3 integration gives a delta function. The x integration 
can then be performed immediately, giving 



L(u) + L(u + tt) = — 2r sin uj cos v , 

1 cosn" 



2r 

M(u) + M(u + tt) = 



k cos v 



(45) 
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It follows that 

t{u) + £(u + 7r) = -2t£; 2 sin u , 
m(it) + m(u + 7r) = 2r(k 2 /k' 2 ) cosu (46) 

and 

+ vr) = x(u) , y(it + 7r) = y(u) + 2rk 2 jk' 2 . (47) 

Thus x, dy/du are strictly periodic functions of it, of period 7r. 

We have derived these periodicity relations for < u < it, but by analytic continu- 
ation they must be true for all real u. 

Evenness / oddness 

We can replace u p in (|l4j) by tt — u p without leaving the domain of validity of the 
equation. Doing so, then negating x, we obtain 

G(tt-u,(3) = G(u,-I3) , (48) 

from which it follows that M(u),m(u),x(u) are unchanged by replacing u by n — u, 
while L(u),£(u),y(u) are negated. Using the above periodicity relations, we obtain 

x(-u) = x{u) , y(-u) + y(u) = -2rk 2 /k' 2 . (49) 

Thus x(u) is an even function, y(u) + rk 2 /k' 2 is odd. 

The above periodicity and evenness/oddness properties agree with those given in 
eqns (5.3), (5.4), (6.9) of Q. 



6 The low-temperature limit k —> 1 



When k — * 1, the branch points (28) pinch onto the real axis and the functions have 
discontinuities in their derivatives at u = (n — 1/2)tt. If —ir/2 < u < vr/2, then 
cost; = cos it and (analytically continuing from < u < vr/2) we obtain 

G(u,0) = {l-e-^+^M/sinh/? , 

A(0) = 2P/(itsmkp) , B{p) = 2/3 2 /(7r 2 sinh/3) . 
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We also find that the constants e, X, Y are 



, N-l 

X = ~2N2 E (N -2j)cot(7r j/N) , Y = 

3=1 



Evaluating L(u), £(u), k' m(u), using the identity 



1 (j-l/2vr-n 
tanu = — V cot '—i , 50) 

N f-; N ' y J 



we find that 

= 2^ E (2i-iv-i) cot ^-V^K - " t y(u) = _ r/ ^ (51) 

3=1 

For ir/2 < u < 2>ir/2 we note that cosv is then — cosu. Some of the terms in the 
above working are then negated and we obtain 



JV+l 

N 1 

x(u 

•> \ - 

3=2 



-L ^(2 J -Ar-3)cot (j ^ " , yH = rA' 2 , (52) 



These results are consistent with the above periodicity and evenness properties. 
Equations (3.25b) and (5.2) of (l| are in error in that the rhs of both should be divided 
by 2. With this correction, (5.2) and (4.13) therein are consistent with ([Sll). 

For — 7r/2 < u < vr/2, we can also write (^) as 

^ sm(7rr/AT) 

from which we can verify that eqn. (5.12) of (lj is correct as written. Fortunately it is 
this equation that is subsequently used in [||, rather than (3.25b) or (5.2), so there is 
no reason to doubt the results of that paper. 
In this limit we can verify that 

A = 2e , (54) 

in agreement with (5.25) of ffl]. 



7 The critical case: k — > 



When k —* the model becomes critical. To investigate the behaviour in this limit 
it is helpful to note that A(/3) and B{(3) can be expressed in terms of hypergeometric 
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functions. Let 

/oo ^2imx ^rj. 

n , , 2 . , 2 x„ ' fc = V COsh • ( 55 ) 
-oo (1 + k z sinn xj n 

Then, taking 9 to be real, we can establish that 

J(m,n) = ~ e 2me (1 + e 2B ) 2n B(n + im, n — im) F(n + im, n; 2n; 1 — e 49 ) , (56) 

provided n is real and positive and |Im(m)| < n. Here B{m,n) = T(m)T(n) /T{m + n) 
is the beta function and F(a, (3; 7; z) is the usual hypergeometric function. We can 
use standard transformation formulae (9.131 of ||) to verify that the rhs of (|56| ) is 
unchanged by negating either m or 9. 

Using formula (9.131.2) of ||, it follows that in the limit when k — ► and 9 — > +00, 

J(m, ti) 5 {e 2ime r(n - im)r(im) + e~ 2imS r(n + im)r(-tm)} /T(n) (57) 

(neglecting terms of relative order e~ or smaller). 

The functions A(j3), B(f3) defined in ( |T8| ) can be written as sums and differences of 
functions J(m, 3/2). Define 

c+(f3) = ^ 3 / 2 r(l + i/3/7r)r(i-i/3/vr) , c_(/3) = n^ 2 T(l - ip/iv)T(± + i/3/vr) . 

(58) 

Then we find that 

A(/3) = k- 1 [e- 2i(ie/7r c + (P) +e 2i/3e/n c.((3)] , 

B(p) = k- 1 \ie- 2il3e l* c + (P) -ie W "c_(/3)l . (59) 



For N = 2, we see from ( |13D that s(f3) = 0, so e, a, A are zero. For N > 2, 
using (|36|), ( p8| ) we can write e as an integral over a quadratic form in A(/3), B((3). 
Substituting the forms (|59|), the explicit /c factors cancel, we obtain terms proportional 
to e _4 * /3e / 7r and e it l 3S / n ; and cross terms where such exponentials have cancelled out. 
To leading order (for k small) the cross terms cancel, so their contribution to e is at 
most of order k 2 . The other two terms differ only by negating (3, so each give the same 
contribution to a and we obtain 



00 



-27t/ s(p)e 4i(3e/7T (3- 1 [c_(/3)] 2 d/3 . (60) 
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We shall also need a: this is of the same form as e, but with an extra factor (2 + 4i/3/-7r) 
in the integrand. 

The integrand in (]60| ) is is an analytic function of (3 on the real axis. The nearest 
singularity in the upper half-plane is a pole at (3 = in/N, coming from s((3), so in 
the limit of 9 large, we can close round the upper half plane and the results will be 
dominated by the contribution from this pole. The residue of s(/3) is (l/87r) sin 2ir/N, 
so we obtain 

e = -rk^ N /{N - 2) , a = -2rk A/N /N , (61) 
where (using standard formulae for the Gamma function) 

r = (2vr 3 )- 1 2- 4 / Ar iV(iV-2)r 2 (l + ^)r 2 (i-^)sin(27r/iV) 

= (iV- 2) tan(vr/iV) r 4 (l/A^)/[4vr 2 Afr 2 (2/iV)] . (62) 

Hence from ([TJ), to leading order (for N > 2), 

A = -A(N-2)rk^~ N ^ N /N 2 , (63) 

in agreement with eqns. (6.2) - (6.4) of HJ. 

From ( |i~7|) and ( p^9|) , neglecting terms of relative order k 2 , 

G(u,0) = sech/3 - ike- l3+2/3u/7T {e- iu - 2i/3e/7T (I + I^/tt)- 1 c + (p) - 

e iu+2if38/* ( i _ 2^/Tr) " 1 c_ (/3) } . (64) 

Now consider the functions L(u),M(u) defined by (p3) as integrals over (3. Substituting 
the above forms of A(f3), B(f3),G(u, /?), we obtain three distinct contributions. 

(i) Terms coming from the sech/3 component of G(u.f3): Using (||), these 
correspond to replacing s(f3)G p (/3) +t(/3) in ( |22] ) by r/(2/3 cosh/3). Their only depen- 
dence on N is via the external factor r, i.e. they are proportional to (N— 1)/N. Using 
(]59|), taking the /3-integral to be principal- valued, and negating f3 in the terms involving 
c+(/3), the integrand contains a factor e 2l/39 ^ T . We can close round the upper half-plane 
and in the limit 9 — > +oo the dominant contributions are from a possible simple pole at 
(3 = and a double pole at (3 = tt/2. The first gives contributions 0, r/k to L(u), M(u), 
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respectively - these cancel the T(k/k' 2 ) cos v term in The second gives contribu- 

tions (di + 2t6/tt) cosn — r(l — 2u/tt) sinu, —(di + 2t9/tt) sinu — r(l — 2u/ir) cos it, 
where <ii is independent of both u and fc. 



(ii) Cross terms: The other contributions come from the terms in (|64| ) involv- 
ing c+(/3), c_(/3). They contain s(/3) as a factor in the integrand, so only occur for 
N > 2. Using (|59D, they can naturally be grouped into three terms: those proportional 
to e- Aif3e / n c + (l3) 2 , c+(/3)c_(/3), and e 4i ^/ 7r c_(/5) 2 . The middle terms are cross terms, 
independent of and k, and give contributions cfecosu, — c^sinu, to L(u), M(u), re- 
spectively, where d2, like di, is independent of both u and k. 

(iii) Terms proportional to fractional powers of k: Finally there are the 
first and third terms mentioned in (ii). We negate the variable (3 in the first term, 
which puts it into a form similar to the third. The resulting integral can be closed 
round the upper half (3 plane. There is no pole at the origin and the nearest pole to 
the real axis is a simple pole of s{(3) at f3 = iir/N. The calculation is now very similar 
to that of e: evaluating the residue we find that the combined contribution of these 
terms to L(u),M(u) is -rk 4 / N (N 2 - 4)" 1 cos u, -rk 4 / N (N 2 - 4)" 1 sinu, respectively. 

Substituting these various contributions into and ( p5|) and working to order 
k 2 log k in terms in x(u) that are independent of u, to order k 2+i l N in other terms, we 
obtain 

ll2 , .,.„. rk 4 l N rk 2+4 / N cos 2u 
x(u) = 2Tr- 1 Tk 2 log(k/2)- w -^ + j^—^ 



. rk 2+4 / N sm2u . . 

y(«) = -T(l-2u/Tr)k 2 j^—^ . (65) 

These results agree with (and slightly extend) eqns. (6.6), (6.8) of 0. 



The scaling region 

In an earlier paper Q we obtained the free energy and the functions A, x(u),y(u) in the 
scaling region near k = 0, considerably extending the above results in a series expansion 
in powers of k 4j/N . We do not with to fully re-derive these results here, but note that 
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if we include not just the contributions discussed above from poles at /3 = iir/N, but 
more generally at j3 = ijn/N, where j = 1, 2, . . ., then 

e = ~Y,k ij/N rj/{N -2j) , a = -(2/N) ^ k^/ N r, 

X = -(4/N 2 )^2(N-2j)k^/ N - 4 rj , (66) 

where 

_ N(N-2j) in (2nj\ T 2( l .i\ F 2(l j\ 

= t^j/N)V\j/N)/T 2 {2j/N) . (67) 

The sums are over positive integer values of j and we should restrict j to be less than 
N/2, since for greater values the contributions are smaller than other terms, of relative 
order k 2 , that we are ignoring. 

Similarly, considering the contributions of such a pole to L(u) , M (u) , x{u) , y{u) , the 
term rk 2+A / N /(N 2 — 4) above is replaced by j rj k 2+i: >' N /(N 2 — 4j 2 ) and we obtain 

x(u) = 2^rk 2 log(fe/2) - jj-^ + E ' J N2 _ 4j 2 

. v-^ j r.j k 2+4j / N sin 2u . . 

y {u) = _ r(l _ 2n/7r)A; 2_^^___ . (6g) 

These results agree with eqn. (B2) of H, and with the relevant H\(k) term in eqns. 
(B5), (B6). 



8 The case N = 2 

When N = 2 the model reduces to the Ising model. As remarked above, the function 
s{(5) then vanishes, causing some simplifications, in particular e = a = A = 0. Sub- 
stituting (^) into (^), performing the (3 integration, then integrating by parts in the 
x-variable, we obtain, for < u < ir, 

_ , , i f°° sinh x dx 

L i u ) = T~ 



47r J-oo sin(u + ix)(l + k 2 sinh 2 x) 1 / 2 



1 1 f°° coshx dx 

^ Akk' 2 ~ 4irk' 2 J-oo sin(u + ix)(l + k 2 sinh 2 x) 1 ' ' 2 ' ^ ' 
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Substituting these into (pi]), (|25|), the factor sm(u+ix) cancels out of the expression 
for x(u), leaving, in agreement with (5.38) of JlJ, 

x( u ) = - k 2 K'/ (2vr cos v) , (70) 

where 

roo 

K' = / (1 + fe 2 sinh 2 x)- 1/2 dx (71) 
Jo 

is the usual complete elliptic integral. 
For y(u) we obtain 

, . k 2 r°° cos u cosh x — ik' 2 sin u sinh x , . 

y( U ) = 2 / — 2~r dx ■ ( 72 ) 

Ank' cosv J-oo sin(-u + ix)^/(l + k 2 sinh x) 
To handle this we introduce Jacobi's elliptic functions of modulus k. Define w, t by 



sinu = snu> , cosn= cnw , sinhx = — isnt , coshx = cut , (73) 



where < w < 2K and t is pure imaginary. Then (72) becomes 



y ^ = / J(w,t)dt , (74) 

8nk' J-iK' 



where 

2 (cn went — k' 2 sn w sn t) 



J(w,t) 



dn w ( sn w cn t + cnwsnt) 
By considering the poles and residues of J{w, t) as a function of t, we can verify that 

( > = (lT + »r) (— ) + ( e + 1) (— ) - 2 e7 (! "> • (75) 

where for instance the first bracketted term on the rhs is the sum of the logarithmic 
derivatives of the Jacobi theta functions H and H\, evaluated with arguments (w+t) /2. 
We can now integrate J{w, t) directly with respect to t (taking care to obtain the 
correct branches of complex logarithms), giving 

It seems that the expression in (5.38) of for y(u) is in error: it should be multiplied 
by 2k and the first occurrence of w p replaced by w p — K. When that is done, it agrees 
with ©. 

We have checked that these N = 2 results have the periodicity and eveness / oddness 
properties discussed above, and are consistent with the k = 1 and k — > limiting cases. 
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9 The inversion and rotation relations 



We have already remarked at the end of section 1 that 

A gp = -A pg , (77) 

and that this is the "inversion" relation || of the chiral Potts model. Let R be the 
automorphism such that ur p = u p + n, vr p = —v p . Then, using ref. 0, 

W q>Rp (n) = W pq (n) , W q , Rp (n) = W pq (-n) , (78) 

so replacing p, q by q, Rp is equivalent to merely rotating the lattice through 90°. This 
leaves the partition function Z is unchanged. From ([LTj) it follows that 

Aq,Rp - A pg = log j qAp - log f pq , (79) 

where 

f pq = {det N W pq {i-j)} 1/N I ^\{W pq {n)^ . (80) 

The aim of this section is to explicitly show that the formula (|IT| ) satisfies (|79|). 

We first need to write f pq (or rather its logarithm) in an appropriate form. It is 
given explicitly in eqn. (3.22) of Q. Using (3.33) therein, we can write it as 

fpq = N^(h pq /kT- 1)/2N 9 vq , (81) 



where 



AT-l 

g pq = Y[{2sm[( Uq -u p + 7Tj)/N]}^ 1 - N y 2N , (82) 

3=0 

1/2 



hpg 



cos[N(9 p - <f> q )/2] cos[N(e q - <f>p)/2] 



\ sm[N(6 q - 9 p )/2] sin[iV(^ - <t> p )/2\ \ 
If u p ,u q are real and < u q — u p < it, then all the Boltzmann weights W pq (n), 
W pq (n) are real and positive, and so are g pq ,h pq . Using the relations (|2|), (|3|), we can 
write h pq as 

sm(N8 p )co S [N(0 q -cP p )/2} 
pg k> sin u p sm[N(6 q -6 p )/2] ' 1 j 

Let us write the variables x p ,y p ,\ p of || as x p ,y p ,\ p . Then, using the definition 
(0) of J p , 

z. _ J4>p f. - J{0 p +k/N) 

dj p — c , y p — c , 



k Xp — 1 k I 'Ap ; k i/p — 1 k \p 5 Jp — \p & ^ , 
A p + A" 1 = (1 + k' 2 + k 2 e 2iu ")/k' , A p - A p 1 = 2(jfe/fc') cos v p 

and 

V = i(l-X p X q )/(X p -X q ) . (84) 
Using the identity (56) of j| , we can establish that 

/oo 
e -P+Wfa^)/*[t(fi) + 8 (p)ca&P/amh 2 0\d0 (85) 
-oo 

provided < Re (u q — u p ) < tt. 

Consider h pq as a function of u p , keeping k and u q fixed. Differentiating and using 
the above relations, we obtain 



d log h pq 



(86) 



du p svci{u q — Up) cos v p 

Provided < u q < tt, it follows that d \ogh pq / du p is analytic in a vertical strip in 
the complex u p plane containing the imaginary axis, and tends exponentially to zero 
as Up — > ±ioo. It can therefore be Fourier analyzed in this strip. Doing this and using 
the definition (0), we find that 

dl J Hpq = - [°° e-P +2 ^-^GJ-(3)d(3 (87) 

du p TT J-oc 

provided — tt/2 < Reu p < tt/2 and < Re (u q — u p ) < tt. 

When Up — > +ioo, then v p — > +ioo, A p — > l/k' and h pq — > iX q e~ 2lVq . When 



v 



-ioo, then v p — > — ioo, X p — > oo and h pq — > —iX q . Integrating ( |S7| ) with respect 



to Up, it follows that 



-P+2f3{u q -u p )/iT 

\ogh pq = -\ogJ q -V / G q (-f3)dp 



oo 



We have introduced several complex variables, but if real and < u a — 

Up < tt, then we recall that G q {—(3) is a real function, and h pq , J q are real and positive: 
(|86|), (|87|), ( |88D are then real equations. The additional conditions —tt/2 < u p < tt/2, 
< u q < tt are needed to ensure convergence of the integrals as written. 

Noting that Xr p = A" 1 , from (^|) and ( |S4"| ) we can verify that 

9q,Rp = 1/ ffpij ; hq,Rp = ^-/h pq , 
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so from (Kl 



log f q ,R p - log f pq = -2 log g pq - 2t log h. 



OO 



= -r\ogJ q + V e~ p+2 ^ u i- u ^^X q (P)d(3 , 

J —oo 

where 

= T (5- l G q {-(5) +2t(fi) + 2s(/3)cosh/3/sinh 2 /3 . (90) 
Using (|43|), we can write this last relation as 

X q (P) = 2s(f3)[G q (-l3) + cosh/?/sinh 2 $ +2t(J3)[l + G 9 (-/?) cosh/3] . (91) 



We can now explicitly verify that the result (11 ) of g, [| satisfies the rotation 
relation (|79|), The relation ( |42|) implies 



G Rp ((3) = -e 2/3 G p (/3) + 2e /3 . (92) 
Using this identity and ([□] ), and noting that Jr p = 1/J p , we find 

-A R -A =- T logJ q +V e-fi+W*,-«p)/* X (p) ; (93) 
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the function G p {[3) having cancelled out of the integrand. The RHS of this result is the 
same as that of (p9|), while from fl77| ) the LHS is Ag 5 ^ p — A p(? : ( |79| ) is therefore satisfied. 

10 Summary 

We have obtained the explicit expressions (|2l| ) - ( |25| ) for x p ,y p . We have shown that 
they satisfy the differential relations (^), and that they agree with previous results 
(corrected where necessary) for the low-temperature case k — * 1, for the critical case 
k — > 0, and for the Ising case N = 2. 

In section 9 we have also explicitly shown that the free energy satisfies the inversion 
and rotation relations. For other models that enjoy the "rapidity difference" property, 
there is a uniformizing substitution that makes it easy to derive the free energy from 
these relations and a simple analyticity assumption. We do not know how to do this for 
the chiral Potts model, nor do we know how to solve similar relations for the generalized 



spontaneous magnetization [l(j, 11, 12]. The presentation in section 9 may cast some 



light on this outstanding problem. 
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